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INTRODUCTION 

It is well known that recent astrophysical data are in agreement with a universe in current 
phase of accelerated expansion, in contrast with the predictions of Einstein gravity in 
FRW space-time. It seems that the most part of energy contents (roughly 75%) in the 
universe is due to mysterious entity with negative pressure: the so dubbed Dark Energy. 
The simplest explanation is ACDM: Einstein gravity plus a small positive cosmological 
constant suffers from the coincidence problem and the cosmological constant issue: 
Kob/\h — 10^^^*^, while if we assume supersymmetry, one has Aob/^th — 10^^^. In 
this case w^e ~ ^ ~ There exist alternative explanations. Among many, we recall: 

i. Modification of gravity on large scale: DGP brane- world model 

ii. Dark energy associated with cosmological scalar fields, quintessence if Wde > 
— 1, phantom matter if Wde < — 1. 

iii. Modified gravity models: R — > R + f{R). 

All these models look like ACDM, but with an effective non constant cosmological 
term. 



MODIFIED GRAVITY AS MODELS FOR DARK ENERGY 



The ACDM model is the simplest possibility but, it is worth investigating more general 
modifications, possible motivations run from quantum corrections to string models: (for 
a recent review see [0] and references therein). We shall consider the modification of 
the kind F{R) = R + f{R). Models of this kind are not new and they have been used in 
the past by many authors, for example as models for inflation, f{R) = aR^ [3]. Recently 
their interest in cosmology was triggered by the model f{R) = —ji^/R, proposed in 
order to describe the current acceleration of the observable universe [0]. 

It is important to stress that these F{R) models are conformally equivalent to Ein- 
stein's gravity, coupled with a self-interacting scalar field, Einstein frame formulation. 
We will consider only the Jordan frame, in which the dynamics of gravity is described 



by F{R) with minimally coupled matter. Observations are typically interpreted in this 
Jordan frame. 

Finally, we would like to mention the so called viable F{R) models, which have re- 
cently been proposed (J], with the aim to describe the current acceleration with a suitable 
choice of F{R) = R + f{R), but also to be compatible with local stringent gravitational 
tests of Einstein gravity = R. The main idea is the so called disappearing of cosmo- 
logical constant for low curvature, and mimicking the ACDM model for high curvature. 
Thus, the requirements are: 

a. f{R) ^ , R^O, compatibility with local tests. 

b. f{R) — > — 2Ao , R +°°, description of current acceleration. 

c. Local stability of the matter. 

As a illustration, we recall a recent example of viable model ^ 

where Rq, and Aq are suitable constants. Its advantages are a better formulation in the 
Einstein frame and a generalization that may also include the inflation era. 

THE DE SITTER STABILITY ISSUE 

The stability of the de Sitter solution, relevant for Dark energy, may be investigated in 
these F{R) models in several ways. We limit ourselves to the following three approaches: 

i. Perturbation of Esq. of motion in the Jordan frame. 

ii. One-loop gravity calculation around de Sitter background. 

iii. Dynamical system approach in FRW space-time. 

We shall briefly discuss these three approaches. We may anticipate that the third one 
can easily be extended to more general modified gravitational models. 

Stability of F{R) model in the Jordan frame 

The starting point is the trace of the equations of motion, which is trivial in Einstein 
gravity R = —K^T, but, for a general F{R) model, reads 

3V^f'{R)-2f{R) +Rf'iR) -R = k^T . 

The new non trivial extra degree of freedom is the Scalaron : 1 +f'{R) = e^^. Requiring 
R—Rq — CST, one has de Sitter existence condition in vacuum 

Ro + 2f{Ro)-Rof'{Ro)=0. 

Perturbing around dS: R = Ro + 5R, with 5R = — ^^/(^^"^ ^X^ one arrives at Scalaron 
perturbation Eq. 

^ ^ 6(l+/(i?o)) 



One may read off the Scalaron effective mass 
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Thus, if m2 > 0, one has stability of the dS solution and the related condition reads 

Rof'iRo) 

If < 0, there is a tachyon and instability. Furthermore, one may show that has to 
be very large in order to pass both the local and the astronomical tests and 1 +f'{R) > 0, 
in order to have a positive effective Newton constant. The same result has been obtained 
within a different more general perturbation approach in Jvl] . 



One-loop F{R) quantum gravity partition function 

Here we present the generalization to the modified gravitational case of the study of 
Fradkin and Tseytlin [l8D, concerning Einstein gravity on dS space. One works in the 
Euclidean path integral formulation, with dS existence condition 2Fo = RqFq, assumed 
to be satisfied. The small fluctuations around this dS instanton may be written as 
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Making use of the standard expansion of the tensor field hjj in irreducible components, 
and making an expansion up to second order in all the fields, one arrives at a very compli- 
cated Lagrangian density L2, not reported here, describing Gaussian fluctuations around 
dS space. As usual, in order to quantise the model described by L2, one has to add 
gauge fixing and ghost contributions. Then, the computation of Euclidean one-loop par- 
tition function reduces to the computations of functional determinants. These functional 
determinants are divergent and may be regularized by the well known zeta-function reg- 
ularization. The evaluation requires a complicated calculation [9] and, neglecting the so 
called multiplicative anomaly, potentially present in zeta-function regularized determi- 
nants (see UM), one arrives at the one-loop effective action, here written in the Landau 
gauge 
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The last term is absent in the Einstein theory. As a result, in the scalar sector, one has 
an effective mass = ^ (^'^~^ —Ro^ ■ Stability requires > 0, in agreement with the 
previous Scalaron analysis, and with the inhomogeneous perturbation analysis [jVO. 



Dynamical system approach 



This approach has been used by many authors [[Ili [1, [12, [H, One works in 
a cosmological setting, namely with a FRW metric, and the main idea consists in 
rewriting the generalized Einstein-Friedman equations in an equivalent system of first 
order differential equations, introducing new dynamical variables Q.i 

ja{t)=vm)). 

Here the evolution parameter has been denoted by t. The critical (or fixed ) points are 
defined by v{Q.q) = 0. The key point is: 

Hartman-Grobman theorem: The orbit structure of a dynamical system in the 
neighbourhood of a hyperbolic fixed point is topologically equivalent to the orbit struc- 
ture of the associated linearized dynamical system, defined by a stability matrix Mq. 

Recall that a hyperbolic fixed point is such that its stability matrix Mq does not have 
vanishing eigenvalues. In other words the theorem states that the flux of a dynamical 
system in a neighbourhood of a hyperbolic fixed point can be continuously deformed 
to the flux of the related linearization. As a result, in order to study the stability of the 
above non linear system of differential Eqs. at critical points, it is sufficient to investigate 
the related linear system of differential Eqs.: 

— dQ(t)=ModD.(t), Mq Jacobian matrix evaluated at t^o 
dt 

The solution of the linearization is well known and the evolution is determined by 
the signs of the eigenvalues of Mq. As a result, the non linear system is stable if all 
eigenvalues of the matrix Mq have negative real parts. 

Let us apply this method to study the stability for F{R) models. Introducing new 
variables are defined by 

o o f(R)-RfiR) o XP 
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the dynamical system equivalent to Einstein-Friedman Eqs. reads 

-^Q.R = 2aR{2-aR)aR-i5 (i-ap-ap), 

= 2aF{2-aR) + {af-aR){i-af-ap), 
-^ap = [2{2-aR)-3{w+\) + i-aF-ap]ap, 

here the evolution parameter is a{t) = \na{t) and w = ^, and the function j8 is /3(i?) = 
^Rf'{R) ■ ^^^t '^^^ ^ complete autonomous system as soon as the quantity /3 can 
be expressed as a function of Q,-. This requires the inversion of ^^(i^^/l^^'* = ji- After 



this inversion, in principle, one has /3 = j5{Q.R,Q.f), and may close the above system. 
The possible problems are: non unique inversions, non trivial domains with divergent 
points, ECT. The non linear algebraic system for critical points is 



= 2aR{2-aR)aR)-i5{\-QF-^p), 

= 2Qf{2-^r) + {^f-^r){'^-^f-^p) 
= [2{2-aR)-3{w+l) + l-^f-Qp]Qp. 

In vacuum p = 0, namely Q.p = 0, and de Sitter critical point existence condition follows 
from the solution Q.r = 2 , Q.f = 1, namely Rq = \2Hq and Rq = Rof'{Ro) — 2f{Ro), in 
agreement with the other methods. In order to investigate the stability of this dS critical 
point (2, 1,0) the associated linear system is 

-^5aR = -4daR+2po5aF + 2po5ap, 

da 

= -2 6Q.R + 5af + 6Q.p, 
da ^ 

—d^p = OdaR + Odap-Sydap, 
da 

and one can read off the stability matrix Mq, whose eigenvalues are Ai = — Sy, 7 > 
and 



A2,3 = ^(-3±v/25-16/3o 



The stability condition associated with the de Sitter critical point requires that the real 
part of all eigenvalues has to be negative, thus 

l+/(^o) . . 
Rof'iRo) ' 

again in agreement with Scalaron perturbation analysis and one-loop de Sitter calcula- 
tion. In the matter-radiation sector, where Q.p is non vanishing, other critical points, in 
general, exist, but their analytical determination, in realistic cases, is problematic, since 
one has to know explicitly /3 in order to close the system, and numerical analysis, in 
general, is necessary. 

We conclude recalling that, within this approach, it is not difficult to deal with 
generalizations of the kind F{R) — > F{R,G,Q, ..) , [15] where F depends on arbitrary 
invariants of tensor curvature as G, Gauss-Bonnet invariant, Q the square of Riemann 
tensor, and so on. In the case G, 2), the associated de Sitter existence solution reads 

F = — ^iFQ + Fq) and the related stability condition is (see [il6i1 and references 
therein) 



R 



^RR^ 3^Q'^ 3^(^GR^^Rq) + l^^^^GG^^^GQ^ ^Qq) 



CONCLUDING REMARKS 



Modified gravity may be seen as the phenomenological description of a fundamental un- 
known theory. From this point of view, corrections to Einstein-Hilbert action depending 
on higher order curvature invariants are likely to be expected (Lovelock gravity is an 
example). 

Among many existing approaches, three methods have been illustrated in order to 
investigate the stability of these models around de Sitter critical points, and the dS 
stability conditions has been derived in all the three approaches. 

These methods have owns advantages and problems, but, in our opinion, the third 
one, the dynamical system approach, permits to study critical points and stability for 
modified gravitational models depending on arbitrary geometric invariants, generalising 
the results obtained for F{R) models. 
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